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; ABSTRACT 

. We have discovered a new rapidly oscillating Ap star among the Kepler Mission target stars, 

KIC 10195926. This star shows two pulsation modes with periods that are amongst the longest 
known for ro Ap stars at 1 7 . 1 min and 18.1 min, indicating that the star is near the terminal age 
main sequence. The principal pulsation mode is an oblique dipole mode that shows a rota- 
tionally split frequency septuplet that provides information on the geometry of the mode. The 
, secondary mode also appears to be a dipole mode with a rotationally split triplet, but we are 

H ■ able to show within the improved oblique pulsator model that these two modes cannot have the 

. . .' same axis of pulsation. This is the first time for any pulsating star that evidence has been found 

for separate pulsation axes for different modes. The two modes are separated in frequency by 
55 //Hz, which we model as the large separation. The star is an CVn spotted magnetic vari- 
able that shows a complex rotational light variation with a period of Piot = 5.68459 d. For 
the first time for any spotted magnetic star of the upper main sequence, we find clear evidence 
of light variation with a period of twice the rotation period; i.e. a subharmonic frequency of 
frot/2. We propose that this and other subharmonics are the first observed manifestation of 
torsional modes in an roAp star. From high resolution spectra we determine Tcff = 7400 K, 
log.g = 3.6 and vsmi ~ 21kms^^. We have found a magnetic pulsation model with fun- 
damental parameters close to these values that reproduces the rotational variations of the two 
obliquely pulsating modes with different pulsation axes. The star shows overabundances of 
the rare earth elements, but these are not as extreme as most other roAp stars. The spectrum 
is variable with rotation, indicating surface abundance patches. 

Key words: Stars: oscillations - stars: magnetic fields - stars: chemically peculiar - stars: 
variables - stars: individual (KIC 10195926) - stars: spots. 
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1 INTRODUCTION 

The rapidly oscillating Ap (roAp) stars are a unique laboratory in 
which the interaction between stellar pulsation and a magnetic field 
can be studied in greater detail than in any other star except the 
Sun. Although the pulsation periods in roAp stars are about the 
same as in the Sun (a few minutes), they are probably driven by the 
K-mechanism in the hydrogen ionization zone, rather than being 
stochastically excited by convection, as in the Sun. It is thought that 
suppression of convection at the magnetic poles reduces damping 
in the hydrogen ionization zone and allows modes of high radial 
order to be driven (Balmforth et al. 2001). 

However, when comparing theoretical instability strips to the 
positions of the known roAp stars in the HR diagram, there remains 
a difficulty in explaining the presence of pulsations in the coolest 
roAp stars. The theoretical red edge seems to be hotter than its ob- 
served counterpart, a fact that was noted when the first theoretical 
instability strip was computed (Cunha 2002), and that is even more 
obvious today, with more roAp stars having reliable effective tem- 
perature determinations below 7000 K. Moreover, the agreement 
does not seem to improve when models with different global metal- 
licity or local metal accumulation are considered (Theado et al. 
2009). While the iron abundance may influence the excitation, it 
seems not to have an impact on the red edge of the theoretical in- 
stability strip. 

The Ap stars have strong, global magnetic fields in the range 
1 — 25 kG and overabundances of some rare earth elements that 
can exceed lO'' of the solar value. These abundance anomalies are 
located in patches on the stellar surface that give rise to photo- 
metrically observable rotational light variations. The abundances 
are also vertically stratified in the observable atmosphere. They are 
thought to arise by atomic diffusion in the presence of a magnetic 
field. Thus Ap stars in general, and roAp stars in particular, are 
test beds for atomic diffusion theory, which is important for so- 
lar physics, stellar cluster age determinations, pulsation driving in 
main sequence and subdwarf B stars and a possible solution for 
frequency anomalies in /3Cep stars. 

The general view of roAp stars is that the pulsation axis is not 
aligned with the rotation axis, thus uniquely allowing the pulsation 
modes to be viewed from vaiying aspect with rotation. This pro- 
vides geometric mode identification information that is not avail- 
able for other asteroseismic targets. A simple model is that the pul- 
sation axis is aligned with the axis of the magnetic field, which is 
assumed to be roughly a dipole inclined with respect to the axis of 
rotation. As the star rotates, the geometry of the pulsation changes, 
leading to amplitude and, in some cases, phase modulation. This is 
the oblique pulsator model (Kurtz 1982). 

In early work it was assumed that the eigenfunction could be 
represented as a pure axisymmetric spherical harmonic, in which 
case the variation of pulsational amplitude with rotation leads di- 
rectly to a determination of the spherical harmonic degree, I. Fur- 
thermore, amplitude variation also allows constraints on the relative 
inclinations of the rotational and magnetic axes to be placed. We 
now know that the eigenfunction is far from a simple spherical har- 
monic (Takata & Shibahashi 1994, 1995; Dziembowski & Goode 
1996; Cunha 2006; Saio & Gautschy 2004; Bigot et al. 2000) and 
values derived in this way should be treated with caution. Treating 
the combined effects of rotation and magnetic field on the pulsa- 
tion modes is complex. Bigot & Dziembowski (2002) showed that 
in the case of a relative weak global magnetic field - Bp ^ 1000 G 
- the pulsation axis is, in general, neither the magnetic axis nor the 
rotation axis, but lies between these. 



There are 40 roAp stars known from ground-based observa- 
tions that are listed by Kurtz et al. (2006), and a few discovered 
more recently (see, e.g., Kochukhov et al. 2009; Elkin et al. 2010). 
High-resolution spectroscopy using large telescopes allows suffi- 
cient time resolution for the detailed study of line profile varia- 
tions caused by rapid pulsation. It turns out that detection of low- 
amplitude roAp pulsations is easier using radial velocities obtained 
in this way rather than high-precision photometry. Several roAp 
stars are known in which the radial velocity variations can be 
detected, but no corresponding light variations can be seen from 
ground-based observations. A global network of spectrographs 
dedicated to asteroseismology, such as the SONG (Stellar Obser- 
vations Network Group) project (e.g., Grundahl et al. 2009) may 
greatly increase the number of known roAp stars. For now, how- 
ever, photometry is still the most cost-effective method for global 
asteroseismic studies in these stars where the frequencies are the 
fundamental data. For detailed examination of pulsation mode be- 
haviour in the observable atmospheres, spectroscopic studies are 
superior (see, e.g., Ryabchikova et al. 2007, Freyhammer et al. 
2009). 

The Kepler Mission has revolutionised our ability to detect 
low-amplitude light variations in stars, reducing the detection limit 
to only a few micromagnitudes even in rather faint stars. The Kepler 
asteroseismic programme is described by Gilliland et al. (2010), 
and details of the mission can be found in Koch et al. (2010). Ke- 
pler is observing about 150000 stars continuously over its mini- 
mum 3.5-y mission lifetime with 30-min time resolution. A small 
allocation (512 stars) can be studied with 1-min time resolution. 
While there were no known roAp stars in the Kepler field be- 
fore launch in 2009, three have been found in the I-min cadence 
data. KIC 8677585 was studied by Balona et al. (2010) who found 
several pulsation frequencies with periods around lOmin, includ- 
ing an unexpected low frequency of 3.142d^^ which is possibly 
a g-mode pulsation. A second roAp star discovered in the Kepler 
data, KIC 10483436, shows a rotational light variation and two fre- 
quency multiplets (Balona et al., in preparation). 

The third roAp star discovered in the Kepler data is 
KIC 10195926 (a = 19 05 27, 5 = +47 15 48, J2000; Kp = 
10.57), which we report on here. It has a principal mode that shows 
a clear frequency septuplet with separations equal to the rotation 
period of the star, Prot = 5.68459 d. Rotational light variations 
are clearly seen, indicating that it is an CVn star. It appears that 
in this star both pulsation poles can be viewed over the rotation 
cycle, which resembles the well-studied roAp star HR3831 (see, 
e.g., Kurtz et al. 1990, Kurtz et al. 1997, Kochukhov et al. 2004, 
Kochukhov 2006, and references therein). A second, much lower 
amplitude mode that shows a frequency triplet split by the rotation 
frequency is also found, thanks to the 5 /imag precision of the Ke- 
pler data. Interestingly, it appears that these two modes do not share 
the same pulsation axis, an effect that has never been detected pre- 
viously. We present theoretical support for this hypothesis. 



2 FUNDAMENTAL PARAMETERS FROM 
SPECTROSCOPY 

Spectroscopic observations of KIC 10195926 were made with 
the high-resolution Fibre-fed Echelle Spectrograph (FIES) at the 
Nordic Optical Telescope (NOT) in 2010 June, July and August. 
The spectra have a resolution of -R = 65 000 and cover the spec- 
tral range 370 — 730 nm. The signal-to-noise ratio (S/N) attained 
is about 60 for the June spectra and 40 for observations in July 
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Figure 1. Top: The long cadence light curve of KIC 10195926 from QO + Ql. The ordinate is BJD - 2450000. Middle: The short cadence light curve of 
KIC 10195926 from Q3.3. Bottom: the short cadence light curve after high pass filtering (described in Section 4) has removed the rotational light variations, 
leaving only the pulsational variations. Only the envelope of the pulsation is visible here, as the pulsation peiiod is too short to be resolved at this scale. Note 
that pulsation maximum coincides with rotational minimum light. 



and August. The star was also observed with the Sandiford Echelle 
spectrograph on the 2. 1-m telescope at McDonald Observatory. The 
spectra have a resolution of i? = 60 000 and were collected on four 
nights in 2010 July covering the spectral region 508 — 605 nm. On 
nights when multiple exposures were taken, they were co-added. 
The S/N for these observations ranges from 50 to 80. Table 1 gives 
a journal of the observations indicating their rotational phase (see 
Section 3 below). 

Most roAp stars have strong magnetic fields in the range of 
1 — 25 kG. Knowledge of the magnetic field strength is impor- 
tant for understanding and modelling these stars. In the spectra of 
KIC 10195926 there is no obvious Zeeman splitting of magneti- 
cally sensitive spectral lines. We therefore searched for magnetic 
broadening by comparing the observations with synthetic spectra 
calculated for the spectral region around the Fe II 6149.258 A line 
using the SYNTHMAG program by Piskunov (1999). With the S/N 
of our spectra and rotational broadening, we can only place an up- 
per limit on the magnetic field modulus of about 5 kG. 

Because of the probable magnetic broadening, it is necessary 
to use lines with negligible Lande factors to measure the projected 
rotational velocity v sin i. Because of the patchy distribution of el- 
ements over the surface, it is necessary to examine spectra at differ- 
ent rotation phases, since lines from elements concentrated in spots 
do not always sample the full rotation velocity of the star From 
the Fe I line at 5434.523 A we determined from spectra at different 
rotation phases that vsini = 21 — 24kms^^. Most spectral lines 
are fitted well with values in this range, but some require an even 
smaller vsini, probably because of the spots. At some rotation 
phases 24kms~^ is needed to fit the line broadening. Our internal 



Table 1. A journal of spectroscopic observations obtained with FIES on 
the Nordic Optical Telescope, and with the Sandiford Spectrograph on the 
McDonald Observatory 2. 1-m telescope. The rotational phase has been cal- 
culated using the ephemeris $rot = 245 5168.91183 -I- 5.68459. 



Nordic Optical Telescope - 


PIES 


JD -2450000.0 


exposure time (s) 


rotational phase 


5369.64536 


1800 


0.31 


5369.66862 


1800 


0.32 


5407.37706 


1800 


0.95 


5407.40009 


1800 


0.95 


5408.37815 


1800 


0.13 


5408.40118 


1800 


0.13 


5411.68126 


1800 


0.71 


5411.70415 


1800 


0.71 


McDonald 2.1 


m telescope - Sandiford spectrograph 


5397.92206 


2400 


0.28 


5398.69240 


1800 


0.42 


5399.70613 


1800 


0.60 


5400.90493 


1800 


0.81 



error on a determination of v sin i is 1 km s^^, which we cautiously 
increase to 2 km s~^. Our estimate is v sin i = 21 ± 2 km s^^, but 
we consider this to be near a lower limit for the equatorial rotation 
velocity, given the larger v sin i that we see at some rotation phases. 

KIC 10195926 has photometrically determined parameters 
from broadband photometry given in the Kepler Input Catalogue 
(KIC): Kp = 10.57, = 7400 K, log 5 = 3.6. The KIC gives 
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an estimate of the contamination factor for this star of 0.015, mean- 
ing that less than 1.5 percent of the light in the observing mask can 
be from contaminating background stars. This makes it exceedingly 
improbable that any of the light variations we discuss in the paper 
might be from a contaminating background star. 

To determine the effective temperature of the star, Balmer 
lines profiles were compared with synthetic profiles starting with 
parameters close to those from Kepler Input Catalogue. Synthetic 
spectra were calculated with the SYNTH program of Piskunov 
(1992). We used model atmospheres from the NEMO database (Vi- 
enna New Model Grid of Stellar Atmospheres) (Heiter et al. 2002). 
This grid has an effective temperature step of 200 K, which we in- 
terpolated to get models with a 100 K step. The spectral line list 
for analysis and SYNTH calculations was taken from the Vienna 
Atomic Line Database (VALD, Kupka et al. 1999), which includes 
lines of rare earth elements from the DREAM database (Biemont 
et al. 1999). 

We found a best fit for Ha and H/3 with synthetic profiles for 
Tes — 7200 K, with estimated error of 200 K. The determination of 
the surface gravity is difficult at this temperature. The Balmer line 
profiles are relatively insensitive to \ogg, constraining it only to 
within ±0.5. lonisation equilibrium for Fe I and Fe II lines and CrI 
and Cr II lines is potentially much more sensitive for determining 
log g, but we did not get reliable results with our initial study. We 
expect a more detailed analysis, which is underway, will better con- 
strain log g. We therefore use here the photometrically determined 
value of log g — 3.6 (cgs) from the KIC. Because the spectrum is 
not very peculiar, we estimate a Icr error of ±0.3. 

A detailed abundance analysis is in progress. We made pre- 
liminary estimates of the abundances of some elements at two 
rotation phases to get an overview of the spectral character of 
KIC 10195926. Using synthetic spectra calculated with SYNTH for 
model Teff — 7200 K and logg — 3.6 we examined some abun- 
dances using the first and third NOT spectra listed in Table 1 . The 
spectral lines are clearly variable with rotation, so other spectra will 
give values that differ for some elements. Surprisingly for a roAp 
star, there are no obvious lines of rare earth elements. The usual 
strong lines in roAp stars of Nd III and Eu II are present and indi- 
cate overabundances compared to the Sun, but less than is typical 
for most other roAp stars. Si II lines are noiinal. There is no obvious 
core-wing anomaly (Cowley et al. 2001) in Ha, a usual signature 
of roAp stars. 

Table 2 gives the abundances we found compared to solar val- 
ues. The second column gives the abundances for a spectrum ob- 
tained at rotational phase 0.31. For this spectrum no clear evidence 
lines of rare earth elements was found. There are some hints of the 
presence of such lines in the spectrum, but the better S/N spectra are 
required. The fourth column gives abundances for a spectrum ob- 
tained at rotational phase 0.95. This spectrum definitely shows lines 
of Eu II and a trace of weak lines of Nd III which confirm the pecu- 
liarity of this star. Lines of strontium are much wider and stronger 
in this spectrum while the barium lines are much weaker when 
compared with the spectrum at phase 0.31. Iron lines are weaker 
and slightly wider with lower abundances. The spectral lines show 
strong variability with rotational period, demonstrating the spot- 
ted surface structure. By comparing with the roAp star HD 99563 
(Freyhammer et al. 2009), we judge that lines of EuII and Ball 
are fomed in rather small spots. Other elements also demonstrate 
nonuniform distribution over stellar surface. Fig. 2 shows the Eu II 
6645-A line at four rotation phases. While better S/N is needed for 
further study, it is clear that Eu II is concentrated is at least two 
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Figure 2. This shows the Eu II 6645 A Hne at four different rotation phases. 
The spectra have been offset on the y-axis for visibility. The bottom spec- 
trum is for rotation phase 0.95, which is close to the time of pulsation maxi- 
mum and also rotational light minimum (see Fig. 4.) This is when one of the 
Eu II spots is close to crossing the line-of-sight. In the middle panels, which 
ai'e near quadrature, it is likely that two small spots of Eu II that are concen- 
trated near the magnetic poles are seen separately. This will be studied in 
detail with future, higher S/N spectra. 

small spots close to the pulsation pole (and hence probably close to 
the magnetic pole) of this star. 

We note that the abundances determined here are not char- 
acteristic of the entire star. Atomic diffusion causes atmospheric 
stratification in Ap stars, with some elements being enhanced by 
orders of magnitude with respect to the Sun, while others are de- 
ficient. As a main sequence A star KIC 10195926 has an age less 
than 1 Gyr, hence is expected to have global abundances somewhat 
higher than those of the much-older Sun. This is a consideration 
when modelling the star. 



3 THE LIGHT CURVES AND ROTATION PERIOD 

KIC 10195926 has been observed by the Kepler mission during 
three "rolls" up to this writing. Following each quarter of a 370-d 
Kepler solar orbit the telescope is rolled to keep its solar panels fac- 
ing the Sun. Choices and changes of targets are made on the basis 
of these "quarters", which sometimes are split into approximately 
1 -month "thirds". Most data (for > 150 000 stars) are obtained in 
long cadence (LC) with integration times of 29.4 min. For 512 tar- 
gets data are obtained in short cadence (SC) mode with sampling 
times just under 1-min. The three rolls for which we have data for 
KIC 10195926 are long cadence in QO (quarter zero), which was a 
commissioning roll of less than lOd, long cadence in QI, and short 
cadence in Q3.3 (25 d during the last third of Q3). Table 3 gives a 
journal of the data. 

Fig. 1 shows the long cadence light curve for QO and QI to- 
gether in the top panel and the short cadence light curve for Q3.3 
in the middle panel. A few outlying points were edited from the 
raw data and a linear trend to correct for instrumental drift was re- 
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Table 2. Chemical abundances for KIC 10195926 and the corresponding 
solar abundances (Asplund et al. 2009). The errors quoted are internal stan- 
dard deviations for the set of lines measured. Columns 2 and 4 give, respec- 
tively, abundances detemiined at two rotation phases, 0.31 and 0.95, using 
Pvot = 5.68459 d (see section 3) . The rotation phase refers to the time 
of pulsation maximum, which coincides with rotational light minimum as 
seen in Fig. 4. Columns 3 and 5 give the number of hnes, Ni, used in each 
case. 



Ion 


log N/Ntot 
4> = 0.31 


Ni 


log N/Ntot 
<t> = 0.95 


Ni 


log N/Ntot 
Sun 


Cl 


-3.6 ±0.2 


2 






-3.57 


Mgl 


-4.4 ±0.1 


3 


-4.3 ±0.1 


4 


-4.47 


Si I 


-4.4 ±0.2 


3 


-4.1 ± 0.1 


4 


-4.40 


Sin 


-4.2 ±0.1 


3 


-4.0 ±0.2 


3 


-4:A9 


Cat 


-5.1 ±0.1 


18 


-5.3 ± 0.2 


10 


-5.66 


Sen 


-8.7 ±0.1 


3 


-9.9 ± 0.2 


2 


-8.85 


Tin 


-6.8 ±0.3 


3 






-7.05 


Cri 


-5.8 ±0.2 


4 


-5.4 ±0.4 


4 


-6.36 


Cm 


-5.7 ±0.1 


6 


-5.7 ±0.1 


8 


-6.36 


Fel 


-4.3 ±0.1 


25 


-4.6 ± 0.1 


29 


-4.50 


Fen 


-4.3 ±0.1 


7 


-4.6 ± 0.1 


7 


-4.50 


Nil 


-5.9 ±0.1 


4 


-6.7 ±0.2 


3 


-5.78 


Srll 


-8.8 ±0.2 


2 


-8.0 ±0.2 


2 


-9.13 


Yn 


-8.9 ±0.1 


2 


-8.8 ± 0.3 


5 


-9.79 


Ban 


-9.0 ±0.1 


3 


-10.4 ±0.3 


2 


-9.82 


Ndlll 






-9.3 ± 0.1 


3 


-10.58 


Eull 






-8.4 ±0.1 


4 


-11.48 



Table 3. A journal of the Kepler data for KIC 10195926. The third column 
is the number of data points after a few obvious outliers have been removed. 
The Q 1 data have had a linear trend removed because of instmmental drift. 
The final column gives the root-mean-square variance in each data set of 
one point with respect to the mean. This includes all variance in noise and 
stellar variability. 



quarter 


cadence 


npts 


time span 
days 


rms variance 
^mag 


QO 


LC 


471 


9.71 


1207 


Ql 


LC 


1621 


33.45 


1190 


Q3.3 


SC 


33 679 


25.07 


1191 



moved from the Ql data. The data are in Kepler magnitude, Kp, 
which is in broad-band white hght. It is obvious from Fig. 1 that 
KIC 10195926 is an CVn star; i.e., it has surface spots that pro- 
duce a variation in brightness with rotation. This is a common fea- 
ture of the magnetic Ap stars and implies that the star has a strong, 
global magnetic field, even though that field is yet to be detected, 
and we were only able to place an upper limit on it of 5 kG in Sec- 
tion 2. In some roAp stars, e.g. HR3831 (Kochukhov et al. 2004) 
and HD 99563 (Freyhammer et al. 2009), the rare earth element sur- 
face spots are concentrated near to the pulsation or magnetic poles. 
These rare earth element spots often correlate with rotational light 
variations, particularly in the blue. With the advent of high preci- 
sion space-based observations we can now see rotational light vari- 
ations from spots in Ap stars that would not have been detected with 
ground-based observations, or at least would have given a poorly 
defined light curve. Such is the case for KIC 10195926 where the 
peak-to-peak range in Kp is only slightly greater than 3 mmag. We 
point out, however, that the rotational light variations of Ap stars 
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Figure 3. Top: This figure plots the standard deviation of one observation 
with respect to a fit of a fifth-order harmonic series at each frequency to 
the data. The procedure breaks down at low frequency, and is meaningless 
at zero frequency, hence the choice of frequency range searched. Bottom: 
This is a higher resolution view of the best-fitting frequency that shows 
there is no alias ambiguity, despite the 158-d gap in the data between Ql 
and Q3.3. 



are strongly wavelength dependent, as the result of flux redistribu- 
tion caused by abundance spots, so that KIC 10195926 may well 
show higher amplitude rotational light variations observed in fil- 
tered light such as Johnson U or B. 

To determine the rotational period of KIC 10195926 we co- 
added the 1-min integrations of Q3 to 30-min integrations and anal- 
ysed all of the data listed in Table 3. There is a 158-d gap between 
the end of the Ql data and the beginning of the Q3.3 data, giving 
rise to possible alias ambiguities in an analysis of the full data set. 
As can be seen from Fig. 1 the rotational light curve is highly nonsi- 
nusoidal, so Fourier analysis is not the best choice for period deter- 
mination. We chose to fit a harmonic series by linear least-squares 
to the data to search for the greatest reduction in the variance in the 
data. 

Fig. 3 shows the results of a fifth-order harmonic series fit to 
the data. There is a clear best frequency at v^ot = 0.17591 d~^ 
with no alias ambiguity. While significant harmonics are present 
up to 15 — 20 in the individual quarters of data, no discernible 
improvement in Fig. 3 can be seen with a higher-order fit. 

We then fitted 20 harmonics of frot = 0.17591 d"^ to the 
full data set by least-squares and then nonlinear least-squares to 
optimise the frequency and make an error estimate. This gave our 
best determination of the rotation frequency of Uyot = 0.175914 ± 
0.000004 d~^, which corresponds to a rotation period of Prot ~ 
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Figure 4. The rotational phase curve for the QO and Ql data for 
KIC 10195926. We have chosen the zero point of the time scale to be 
to = BJD 2455168.91183, which is the time of pulsation maximum for 
the principal mode in this star. The diagram shows that pulsation maximum 
coincides with an extremum of the rotational light curve. 



5.68459 ± 0.00013 d. This error estimate is internal. Some system- 
atic effects are still present in the data, such as small instrumental 
drifts and slight changes in amplitude between quarters. These will 
be improved in future data reductions using more detailed pixel in- 
formation. Since more data for KIC 10195926 are being obtained, 
the rotation period will be determined to higher accuracy in future 
studies. The precision given here is excellent for comparison with 
the pulsation phases within the oblique pulsator model. 

Fig. 4 shows the long cadence light curve for QO and Ql 
phased with the rotational period. The double-wave character of the 
rotational light curve indicates two principal spots. The zero point 
of the time scale has been selected to coincide with the time of pul- 
sation maximum, as determined in Section 4 below. It can be seen 
that minimum rotational brightness coincides with pulsation maxi- 
mum, as is expected in the oblique pulsator model (Kurtz 1982), in 
the case when the spots are associated with the magnetic/pulsation 
axis of the star. This coincidence can also be seen by comparing the 
middle and bottom panels of Fig. 1. 

The asymmetry of the light curve requires differences in spots 
on opposite hemispheres. It does not appear that the spots seen at 
rotational phase 0.5 are coincident with the second pulsation pole. 
The spots of Ap stars are often much more pronounced than this 
in filtered light (say, UBVRI), but are generally not studied pre- 
viously in white light. An exception is a Cir, which was studied in 
white light with the star-tracker on the WIRE mission (Bruntt et al. 
2009). While the rotational light curve of a Cir is more symmetri- 
cal, it has a similar low amplitude of only a few mmag. 

Interestingly, Bruntt et al. (2009) showed that the pulsation 
pole for aCir is not coincident with the spots. As that is the case 
here also for the second spot at rotation phase 0.5 of KIC I0I95926, 
we speculate that as more precise rotational light curves are ob- 
tained for Ap stars, we will find more low amplitude rotational vari- 
ation where the spots are not coincident with the magnetic poles. 
We also speculate that this will show a correlation with magnetic 
fields that are less strongly dipolar than in Ap stars with higher am- 
plitude rotational light variations. 

3.1 The rotational subharmonic 

When a 20*'^-order harmonic series is fitted to the rotational light 
curves for the individual quarters, both QO and Q3.3 show a highly 
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Figure 5. Light variation after removal of rotational modulation. The 
time scale is the same in both plots with time zero chosen to be 
BJD 2454950.000. The top panel shows the short-cadence data. Note the 
higher scatter at light minima and light maxima due to the roAp pulsation. 
The bottom panel shows the long-cadence data after correction of long-term 
drift. In both cases a clear variation with twice the rotational period can be 
seen. 



significant peak at ^u^ot, the subharmonic. That this appears in two 
independent data sets shows that it is not an artefact. The residual 
light variation after removal of the 20* ^ -order harmonic fit to the 
rotational light curve is shown in Fig. 5 where it is evident that there 
is periodic variability on a time scale of 2Pi.ot ~ 11.3692 d. It 
is this that gives rise to the subharmonic and its harmonic in the 
frequency spectrum. 

Ellipsoidal variables often have frequency spectra where the 
highest peak occurs at 2i/'rot as a consequence of a double-wave 
light curve with different maxima and minima. If the highest peak 
in the frequency spectrum were mistaken for Viot, then there would 
appear to be a subharmonic present. This is most unlikely to be the 
case for KIC 10195926. The light curve seen in Fig. 1 bears little 
resemblance to an ellipsoidal variable, and the incidence of short 
period binaries in magnetic Ap stars is very low. We are gathering 
high resolution spectra of this star to study its rotational variations. 
Those will be used to test for binary radial velocity variations. It is 
clear that the dominant rotational light variation is of the CVn 
type, i.e. that the star is a spotted magnetic variable. If the variation 
with 2Prot were to be interpreted as orbital, then the rotation period 
would have to be half that of the orbital period, which is again 
unlikely. 

Subharmonics are seen in stellar pulsation, most strikingly in 
the Kepler data for RRLyr stars, including RRLyrae itself (Szabo 
et al. 2010; Kolenberg et al. 2010). This is interpreted in terms of 
nonlinear effects in the pulsation, but this explanation does not ap- 
ply here because we are dealing with rotation, not pulsation. 

We interpret the rotational light variations in Fig. 1 to be 
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Figure 6. The upper panel is an amplitude spectram of the Q3.3 short ca- 
dence data showing the low frequency spot variations. While it cannot be 
seen easily at this scale, the highest noise peaks have amplitudes less than 
5 /xmag. A high pass filter was thus ran to remove all variations at frequen- 
cies less than 0.05 mHz with amphtudes greater than 5 /^mag. The bottom 
panel shows the amplitude spectram in the same range after the high pass 
filtering. Note the change from mmag to fimag between these panels. 
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Figure 7. Top: An amplitude spectram over the full range of pulsation fre- 
quencies. The principal frequency multiplet ai'ound v\ and its haiTnonic can 
be seen, as well as the frequency of a second mode close to the principal at 
V2 ■ Note that the ordinate scale is in /imag. Bottom: An amplitude spectram 
over the full range of pulsation frequencies for the residuals to the fit in Ta- 
ble4. One significant frequency is present at 1/3 = 0.365 mHz (31.54d^^) 
with an amplitude of 6.5 ± 1.2 /^mag. This is formally significant at the 
5.4(7 level, but it is a known Kepler artefact. 



caused by surface spots and we conclude that the rotation period of 
the star is Prot = 5.68459 d. It is unlikely for an CVn spotted 
rotator to have two sets of spots that are as asymmetric as those of 
KIC 10195926, but nearly identical on two opposite hemispheres, 
as would be required if the rotation period were twice the value 
given above. There is no known case of this amongst the CVn 
variables. 

We fitted by least-squares a harmonic series for ^v^ot- Other 
than i^rot and ^Vrot, only the harmonics of v^ot were significant. 
Thus if the subharmonic were the rotation frequency, we would 
have a harmonic series describing the data that included only every 
other harmonic. This, too, is unprecedented and unlikely. 



4 THE PULSATION FREQUENCIES 

Fig. 6 shows an amplitude spectrum of the short cadence data seen 
in the middle panel of Fig. 1. To study the pulsation frequencies 
we ran a high pass filter to remove completely the rotational light 
variations, and any instrumental drift. This was done by sequen- 
tial automatic prewhitening of sinusoids until the noise level at low 
frequencies matches that at higher frequencies. The highest ampli- 
tude noise peaks are less than 5 /imag, so all peaks with amplitudes 
greater than this at frequencies less than 0.05 mHz were removed 
from the data. The purpose of this high pass filter is to ensure white 



noise so that the error estimates for the high frequency peaks are re- 
alistic. There is no significant crosstalk between the low frequency 
rotational variations and the high frequency pulsational variations. 

The bottom panel of Fig. 1 shows these high-pass filtered data, 
where it can be seen that the pulsation amplitude is modulated twice 
per rotation period. This is because the principal pulsation mode is 
a distorted dipole mode and the rotational inclination, i, and mag- 
netic obliquity, /3, are such that i + P > 90°, so that both poles are 
seen. That the maximum amplitude is different when alternating 
poles are presented suggests that neither i nor /3 is close to 90°. 

Fig. 7 shows an amplitude spectrum over the full frequency 
range of variations seen in the star. The principal pulsation mode 
is seen at vi — 972.58 /iHz, and its harmonic can be seen at 
twice that frequency. Both of these are multiplets split by the ro- 
tation frequency. There is another pulsation mode close to ui at 
1^2 = 919.55 /iHz. These are the second lowest frequencies (i.e., 
the second longest pulsation periods) detected in roAp stars, after 
HD 116114 (Elkin et al. 2005), as a consequence of a relatively ad- 
vanced evolutionary state of KIC 10195926. We return to this point 
in Sections 4.6 and 4.7 (see Fig. 16). We also note that the bright 
roAp star /3 CrB (Kurtz et al. 2007) has a similar low pulsation fre- 
quency and also shows overabundances of Nd III and Pr III of only 
1 dex as opposed to the typical 3 dex seen in many roAp stars - 
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Figure 8. Amplitude spectra in tlie frequency range of v\ and V2 ■ Tlie top 
panel shows the multiplets of v\ and V2 ■ The middle and bottom panels 
show them individually at higher resolution. Note the difference in the mul- 
tiplet amplitude patterns, as this is significant for mode identification. 



one of their characteristic signatures (Ryabchikova et al. 2004). We 
return to this point in Section 5. 

Fig. 8 shows a higher resolution look at Pi and U2. The middle 
panel shows the multiplet of v\ . It is a frequency septuplet split by 
exactly the rotation frequency. This pattern is reminiscent of the 
well-studied roAp star HR 3831 (see, e.g., Kurtz et al. 1997; Kurtz 
et al. 1990). It is an indication of a perturbed dipole mode for which 
i + /? > 90°, consistent with the rotational variation seen in Fig. 1 
where there is a double wave. The bottom panel of Fig. 8 shows 
the frequency triplet of P2- This is different to the pattern of v\, as 
here the central peak has the highest amplitude. Thus p\ and V2 are 
either associated with modes of different degree, /, or have different 
pulsation axes. We return to this point in Section 4.3. 
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Figure 9. Schematic amplitude spectrum for the u\ septuplet and ^2 triplet 
showing the component frequencies from the bottom two panels of Fig. 8. 
The splitting between the components of the multiplets is equal to the ro- 
tation frequency. The actual pulsation frequencies are those for the central 
peaks of the multiplets: v\ = 974.6183 /^Hz and U2 = 919.5464 /^Hz. 



Table 4 lists the values of the derived frequencies from a non- 
linear least-squares fit to the data. We found the highest peak in the 
amplitude spectrum, then fitted it, along with all previously identi- 
fied peaks, by linear least-squares to the data. We then prewhitened 
the data by that solution and searched the residuals for the next 
highest peak until the noise level was reached. After all of the fre- 
quencies shown in Table 4 were found, they were then fitted by 
non-linear least-squares to the high-pass filtered data set. As the 
noise is white for these data, the formal errors are appropriate. Af- 
ter prewhitening by the solution in Table 4, there is one significant 
peak left at 31.5 d~^ - a known artefact in the Kepler data - as 
is shown in the bottom panel of Fig. 7. Otherwise there is nothing 
but white noise with highest peaks of 6 /imag. Fig. 9 shows the fre- 
quency multiplets for u\ and U2 schematically. 



4.1 The principal pulsation mode and the oblique pulsator 
model 

The separation of all the multiplets is consistent with the rota- 
tion frequency determined from the rotational light variations, as 
is expected in the oblique pulsator model. We therefore adopt 
the rotational frequency determined in Section 3, which is Prot = 
2.03604 ± 0.00005 fMz. We tested the hypothesis that the rotation 
period may be 2Piot, by fitting by least-squares a 13-multiplet split 
by the subharmonic frequency, ^Vrot, as the v\ multiplet. None 
of the additional peaks between those seen in Fig. 9 has an ampli- 
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Table 4. A non-linear least-squares fit of the ui septuplet and its quintuplet harmonic, plus the U2 triplet. The last column shows the frequency difference 
between the frequency on that line and the one on the previous line. This is the rotation frequency. The value of Uiot = 2.034 ± 0.005 fiHz translates to 
^rot = 5.690 ± 0.014 d, which is consistent with Pi-ot = 5.68459 ± 0.00013 d derived in Section 3 from the rotational light variations. The zero point for 
the phases is BJD 245 5169.0. The residuals per observation for this fit equal 169.8 fimag. 
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.992 ± 0.094 



V\ - V2 



55.072 ±0.028 



tude more than 2.5 /imag, hence all lie well below the highest noise 
peaks of 6 /imag. There is thus nothing in the pulsation modulation 
with rotation to suggest that 2Prot might be the rotation period. 

The frequency septuplet of v\ of KIC 10195926 is very similar 
to that of HR 3831 (Kurtz et al. 1990, 1997), which is interpreted 
as arising from a distorted dipole mode. Kurtz (1992) showed that 
the mode in HR3831 is primarily an oblique dipole mode from 
a spherical harmonic decomposition of its frequency septuplet. 
Kochukhov (2004) found that both dipole and octupole components 
are needed to describe the radial velocity variations in HR3831. 
These results do not imply that both dipole and octupole modes are 
excited. There is only one mode, but the magnetic distortion of that 
mode from a normal mode can be described with multiple spherical 
harmonic components. Dziembowski & Goode (1996), Bigot et al. 
(2000), Cunha & Gough (2000) and Saio (2005) show theoretically 
how the effect of a dipole magnetic field leads to this distortion 
of the pulsation mode. Saio (2005 - his Fig. 12) specifically shows 
a comparison of the latitudinal and horizontal displacements from 
his model with the observational results of Kochukhov (2004) for 
HR3831. 

However, as in the well known roAp stars HR3831, HR 1217 
and aCir, the frequencies of the modes in KIC 10195926 appear 
dominantly to be triplets in photometry. Whether the latitudinal 
magnetic distortion is weak or strong, the mode geometry can be 
reasonably described by a dipole mode, since the higher order com- 
ponents (such as ? = 3) have low visibility when averaged over the 
disk. For a dipole triplet, the phases of all three frequencies should 
be equal at the time of pulsation maximum for a simple oblique 
dipole mode, i.e. <^_i = (j!>+i = <^o at that time. In Table 5 we 
have chosen U = BJD 245 5168.91183 to force = 0+1. 
We then note that differs from the value of and by 
only 0.22 ± 0.02 rad. While this is significant (and is a result of 
the magnetic distortion of the mode), it is close enough to sup- 
port an approximation of the central triplet as arising from a sim- 
ple oblique dipole mode. We also note that Fig. 4 shows that this 



time corresponds to the rotational light minimum, showing that the 
abundance spot that gives rise to this minimum is in the plane de- 
fined by the rotation and pulsation poles. This is confirmed in the 
phase curve of the principal pulsation mode v\ shown in Fig. 10. 
We found the amplitude and phase of this mode by fitting v\ to 
consecutive 5 cycles of data. We have plotted that against rotation 
phase using the same zero point as in Fig. 4. The amplitude modu- 
lation and TT-radian phase reversal are characteristic of an oblique 
dipole pulsation mode (Kurtz 1982). 

The magnetic field present in roAp stars partially raises the 
{21 + 1) degeneracy of non-magnetic modes whose properties (fre- 
quency, amplitude, phase) then depend on the absolute value of 
their azimuthal order m. The modes can be either aligned or per- 
pendicular to the field axis. However, in roAp stars, rotation is 
enough to play an important role in the mode properties. Indeed, 
since the magnetic and rotation axes are in general not aligned 
(/? 7^ 0), the rotation acts as an asymmetric perturbation on the 
magnetic modes, which are thereby coupled with respect to the az- 
imuthal order (rotation is not strong enough in roAp stars to couple 
modes of different degrees /)• 

In this case of rotationally coupled magnetic modes the dipole 
mode displacement vector describes an ellipse during the pulsation 
period which is defined by its two axes: the X-axis in the plane de- 
fined by the rotation and magnetic axes and the other, the Y-axis, 
perpendicular to this plane. This is the generalized oblique pulsator 
model presented by Bigot & Dziembowski (2002). The mode is 
then fully described by the inclination of the ellipse that we defined 
here by the inclination 7 of the X-axis with respect to the rotation 
axis and the angle ?/) whose tangent is the ratio between the lengths 
of the Y- and X- axes; see Bigot & Dziembowski (2002) and Bigot 
& Kurtz (2010) for more details. For a pure m = dipole mode 
in the magnetic reference system we have 7 = /? and ?/j = 0. For 
a pure m = ±1 mode in the same frame we have 7 = 7r/2 + /3 
and -i/i = ±7r/4. The light curves (or the triplet in the frequency 
domain) are described by only three angles: i, 7, V' (Bigot & Kurtz 
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Table 5. A linear least-squares fit of the ui septuplet and its quintuplet 
harmonic, plus the V2 triplet. The rotation frequency has been forced to 
be equal to 2.034 /xHz for the splittings of the multiplets. This is to test the 
phase relationships. The time zero point is BJD 245 5 168.9 11 83; it has been 
chosen to force the phases (/>(vi — frot) = (/'(I'l +frot ), Notably, the phase 
(/((vi) is close to the same value. Exact equality and only a frequency triplet 
is what is expected for a pure oblique dipole with no magnetic perturbation. 
The residuals per observation for this fit equal 169.8 /imag. 



2010). The last two angles are functions of the magnetic field, ro- 
tation, P and frequency. 

The strength of the centrifugal force coupling of the magnetic 
modes depends on the frequency difference between the magnetic 
modes m = and m = ±1 and the centrifugal shift of frequency. 
In order to measure the effects. Bigot & Dziembowski (2002) de- 
fined the parameter to be the ratio between the magnetic 
frequency shift divided by the centrifugal shift: < 1 corre- 
sponds to a dominating rotation regime, = 1 corresponds 
to strong coupling, > 1 corresponds to a magnetic regime 
where there is almost no coupling. The Coriofis force is far weaker 
than the centrifugal force in roAp stars and almost does not enter 
into the balance that defines the inclination of the mode. However, 
it plays a role in the asymmetry of the amplitudes A+\ and A-\. 

We do not know the magnetic field (Bp, P) configuration for 
this star; thus far we have only been able to limit the polar field 
strength to be less than 5 kG. Because rotation is not strong enough 
to couple modes of different I, the presence of latitudinal distortion 
of the mode, seen as a septuplet for i^i , is only due to the magnetic 
field and can be used as a constraint on the field strength. This is 
what is done in the following Section 4.2. The effect of rotation is 
taken into account in Section 4.3. 



4.2 A pure magnetic model for the principal pulsation mode 

Using the formalism of Saio (2005) we have solved for the effect 
of a dipole magnetic field on the i/i mode. In this formalism the 
mode is assumed to be a pure m — mode with its axis of sym- 
metry aligned with the magnetic axis; i.e. 7 = /3. We have run a 
set of models with appropriate Tefi, logg. Unperturbed evolution- 
ary models were calculated by the OPAL opacity with a standard 
chemical composition of {X, Z) = (0.7, 0.02). Convective energy 
flux in the envelope is neglected assuming a strong magnetic field 
suppresses the convection. Helium is assumed to be depleted above 
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Figure 10. This shows the amplitude and phase modulation as a function 
of rotational phase for i/i . It has been created by fitting 5 cycles of data by 
linear least squares for v\. Notice the tt radian reversal at quadrature; this 
is a typical sign of dipole pulsation in an oblique pulsator The difference 
in amplitude of the two maxima is a function of the different aspects as the 
pulsation pole crosses the stellar meridian. Rotation phase zero is the same 
as in Fig. 4, .i.e., with respect to to = BJD 2455168.91183. 



the first ionization zone, treating it in the same way as in Saio et al. 
(2010). 

Here, disregarding the effect of rotation, we assume the pulsa- 
tions to be axisymmetric with respect to the axis of a dipole mag- 
netic field which is inclined to the rotation axis by an angle of /3. 
The eigenfunction of a pulsation mode is expressed by a sum of 
terms proportional to Legendre functions with different values of I. 
Solving numerically the differential equations for linear non- adia- 
batic oscillations including the effects of Lorentz forces, we obtain 
an eigenfrequency and the corresponding eigenfunction; i.e., the 
relative amplitude of each component of the expansion as a func- 
tion of the depth in the envelope, which determines the latitudinal 
dependence of the pulsation. The latitudinal dependence of the lu- 
minosity perturbation on the surface predicts the rotational modu- 
lation of the pulsation amplitude (see Saio & Gautschy (2004) and 
Saio (2005) for details). 

Good fits to the observed amplitude structure seen in Fig. 9 
and the amplitude and phase variations as a function of rotation 
seen in Fig. 10 were found with a model with M = 1.70 Mq, 
logL/Le = 0.998, logTcff = 3.854, logiJ/Ro = 0.314 and 
log g = 4.04. Two good matches to the observations were found, 
one with Bp = 0.7 ± 0.2 kG and the other with Bp = 15 ± 2 kG. 
Our high resolution spectra clearly rule out a polar field strength 
as large as 15 kG, since that would produce clear Zeeman splitting. 
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Table 6. We list all possible values of (i, /3) that lead to an inclination of the 
dipole mode v\ in agreement with the amplitude ratio defined by the data in 
Tables, i.e. (A+i + A_i)/j4o = 3.96 it 0.07. We assume equal magnetic 
and centiifugal strength = 1. Columns 1 and 2 Hst the values of i 

and /3, respectively. Column 3 lists the value of the inclination 7. Column 
4 lists the inclination of the dipole axis with respect to the magnetic field. 
Column 5 lists the corresponding polarization (i/i = is linear). Columns 6 
and 7 list the angles made by the pulsation axis with the line of sight at the 
time of amplitude maxima. Column 8 lists the quantity "D defined in Bigot 
& Dziembowski (2002) to measure the alignment with the magnetic field 
axis (I? = means alignment). The rows in bold characters corresponding 
toi = [62°, 65°] lead to an amplitude ratio (A+i - A_i)/(A+i + 
equal to the observed one. All angles are given in degrees. 
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which is not observed, leaving only the model with Bp = 0.7 kG. 
For the latter model, agreement with the observed pulsation ampli- 
tude and phase variations requires that 7 = /3 = 66° and i — 62° . 

Fig. 1 1 shows the model amplitude and phase variations as a 
function of rotation phase for an n = 16 dipole mode, in good 
agreement with the observations shown in Fig. 10. Fig. 12 compares 
the observed amplitudes of the multiplet structure with the observa- 
tions. A higher than normal value of the limb-darkening coefficient 
of /X = 0.8 was needed to get good agreement for the A±2 and 
^±3 components. The eigenfunction shown in the bottom panel is 
significantly deformed from the Legendre function Pi (cos 8), even 
though A±2 and A±3 are much smaller than Ai and ^o- This is 
because for an odd mode such as a dipole mode, I ^ 3 compo- 
nents (which produce the A±2 and A±3 sidelobes) are visible only 




= 1, n = IB, Bp = 0.700, i/= 0,9789 mHz 
li= 0.80, 7= 66.0, i = 62.0 
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Figure 11. This shows the model pulsation amplitude and phase as a func- 
tion of rotation phase for an n = 16 dipole mode with a frequency of 
u = 978.9 fiHz, close to the observed ui = 974.6 fiHz, and with i = 62°, 
7 = 66°, in good agreement with the oblique pulsator model constraints 
given in Table 6. 



M=1.70, log L=0.996, log T^,j=3.B54, log g=4.Q40, Bp = 0.7 kG 




Figure 12. Top: A schematic amplitude spectrum for the model frequen- 
cies for 1/1. The solid lines show the model amplitudes, while the dotted 
lines (displaced in frequency for visibility) show the observed amplitudes. 
Bottom: The latitudinal dependence of the eigenfunction. The dotted line 
shows pure dipole eigenfunction. 



through the limb-darkening effect; without limb- darkening these 
components would not be visible due to cancellation on the surface. 
This increased limb-darkening coefficient is consistent with the at- 
mospheric temperature gradient in roAp stars, which is steeper than 
in normal stars, as is usually shown by the core-wing anomaly in 
the Balmer lines (Cowley et al. 2001) for most roAp stars. 
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4.3 Rotational distortion of tlie magnetic mode vi 

The magnetic field derived in the previous section is rather weak 
for roAp stars and its effect on the mode inclination can be bal- 
anced by the centrifugal force. In this case, as shown by Bigot & 
Dziembowski (2002), there is no reason to have an axis of pulsa- 
tion aligned along the magnetic axis or the rotation axis. In order to 
calculate the mode inclination 7 and its polarization V', we assume 
a value of = 1 (ratio of magnetic to centrifugal shift of fre- 

quency). The ratio x between the Coriolis and centrifugal shift is 
proportional to ~ n~'^f7~^, where n is the number of radial nodes 
of the mode. Since v\ has a lower frequency and lower n than in 
HR 383 1 , the value of x is larger than for HR 383 1 . From the stellar 
model and eigenfunctions used in Section 4.2, we find x = 0.045. 

Determining the mode inclination from light curves depends 
strongly on the viewing aspect of the mode, i.e. the value of i. In- 
deed, we found several solutions (i, /?) leading to a fit of the ampli- 
tude ratio (^+1 + A_i)/Ao, as long as i < 72° - see Table 6. The 
values follow roughly the law i + /3 « 80 — 100° . While the angles 
derived in Table 6 lead to the same goodness of fit to the light curve, 
they lead to completely different values of the dipole inclination 7. 
Indeed, low values of i lead to a dipole almost aligned with the 
magnetic axis, whereas large values lead to a very inclined mode: 
e.g. 7(i = 19°, /3 = 80°) = 85°, 7(1 = 41°, /3 = 65°) = 77° 
and 7(1 = 62°, /3 = 37°) = 64°. In all cases, the polarization 
remains small showing that this mode is almost linearly polarized. 
Due to the uncertainty in i, it is hard to reach a conclusion about 
the inclination of the principal mode in KIC 10195926. The same 
problem appeared in the case of HR3831. Bigot & Dziembowski 
(2002) found a mode well inclined with respect to the magnetic axis 
7 — /3 « 26° based on a value of i = 89° available from spectropo- 
larimetry at this time. A more recent and secure value of i = 68° in 
HR3831 led Bigot & Kurtz (2010) to a mode almost aligned with 
the magnetic axis 7 — /3 ~ 0. The alignment between magnetic 
and pulsation axes in HR3831 seems natural regarding the large 
magnetic field found for the star ~ 2.5 kG. In KIC 10195926 the 
lower value of Bp — 0.7kG derived in Section 4.2 suggests a more 
inclined mode to the magnetic field. 

We can extract more constraints on the inclination by consid- 
ering the inequality of the side peaks ^+1 7^ A _i. The origin of the 
inequality is due to the Coriolis force which, unlike the magnetic 
and centrifugal forces, does not act in the same way on the ±m 
components of the modes. The ratio {A+i — A-i)/{A+i + A-i) 
is therefore very sensitive to the value of x which determines the 
polarization of the mode (Bigot & Dziembowski 2002). It also de- 
pends on the magnetic field and centrifugal forces. Using the con- 
straint on the unequal amplitude ratio, we restrict the allowed val- 
ues of the rotational inclination and magnetic obliquity to i ~ 63° 
and P ~ 34° . The mode is well inclined to the magnetic axis with 
7 — /? ~ 26°. We used the formalism developed in Bigot & Kurtz 
(2010) to fit of the light curve of ui. An example is shown in Fig. 13 
for(i = 62°,/3 = 63°). 

4.4 The observed amplitude 

The peak pulsation amplitude of KIC 10195926 occurs at rotation 
phase zero in Fig. 1. The highest amplitude occurs when the cen- 
tral three frequencies of the septuplet are in phase, which gives 
{A^:}1 + A'o^'^ + A'^^1) = (168 + 79 + 144) ^mag ?a 0.4mmag, 
as is seen in Fig. 10. This is potentially detectable in ground-based 
observations, but not easily. However, this is the Kepler white light 
bandpass amplitude. A study of the roAp star a Cir using the WIRE 
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Figure 13. Amplitude (top panel) and phase (bottom panel) of the light 
curve variation coiTesponding to the mode u\ as function of the rotation 
phase. The full line corresponds to a mode tilted to the magnetic axis by 
27° with a small polarization (-4°); i = 63°. 

satellite star-tracker by Bruntt et al. (2009) in comparison with si- 
multaneous ground-based photometry through a Johnson B filter 
found a weighted mean amplitude ratio of j4s/Awire = 2.3±0.1. 
If we assume that the ratio Ab / AKepier is similar, then the maxi- 
mum amplitude that would be observed for KIC 10195926 in B is 
expected to be 0.9 mmag, or peak-to-peak nearly 2 mmag. There 
are many other roAp stars known with similar small amplitudes 
(see Table 1 of Kurtz et al. 2006). 

Columns 6 and 7 of Table 6 give the angle between the line- 
of- sight and the pulsation pole for various possible geometries for 
KIC 10195926. These show that an implausibly high intrinsic am- 
plitude greater than 20 mmag - more than is observed in any roAp 
star - would be needed if i were less than « 10°, so those geome- 
tries can be ruled out since 1/ cos(a(0)) and 1/ cos(a(7r)) are too 
large (> 10). 

From the spectroscopically derived values we found in Sec- 
tion 2 of Tcff = 7200 K and log 3 = 3.6, we estimate that 
R — 3.6 Rq. With the known rotation period Prot = 5.6836 d, 
this radius gives v^ot = 32kms^^. With the measured vsini — 
21kms^^, we then find i ~ 40°. Although there are substantial 
uncertainties in this estimate, the angles in Table 6 suggest that we 
are seeing nearly the full intrinsic amplitude of KIC 10195926 at 
rotation phase of best aspect. 

If we use the radius and inclination of our best-fitting model 
from Section 4.2, R — 2A Rq and i — 62°, then we find that we 
expect u sin i — 17 km s ^ ^ . While the errors involved are relatively 
large, this suggests that the model radius and gravity are too low, 
and the spectroscopically derived values from Section 2 are more 
probable. 

4.5 Tlie harmonic frequencies of ui 

Fig. 14 shows a schematic amplitude spectrum for the quintuplet 
centred on 21/1. Since the frequency of the central component is 
twice the frequency of the central component of the principal mode 
septuplet, it appears to be the first harmonic of the principal pulsa- 
tion mode vi. Since the principal mode vi is also linearly polarized 
along a fixed axis in the plane defined by the magnetic and rota- 
tion axes, we can consider that it is an m = mode with an axis 
inclined by 7 with respect to rotation axis. 
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Figure 14. Schematic amplitude spectrum for tlie 2v\ liamionic quintuplet. 



We therefore use the formalism developed by Kurtz et al. 
(1990) for nonlinear dipole oscillations in roAp stars and replace 
P by 7. We then expect 
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in Table 4 we find for 
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1, which is consistent with the observations. 
The amplitude ratios between the central component of the 
quintuplet of 2v\ and the side components are given by Eqs25 
and 26 of Kurtz et al. (1990). If we define by 2; = tan i tan 7 



then we have 
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(2fl) 
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(2-'l) 



/A\,^' = x^/{x^ + 2) and 

1+2 + ^-2 ' ' ) Mc?"' ' = 4a-/ (a-^ + 2) . We find for both ra- 
tios ~ 0.89, which is again fully consistent with the observations. 

It is reasonable to conclude from these considerations that the 
quintuplet centred on 2u\ is indeed the first harmonic of the princi- 
pal pulsation mode vi . It is potentially possible to derive the values 
of i, (5 and /i from the amplitude ratios of the multiplets of the 
principal mode v\, the secondary mode V2, and the harmonic 2v\. 



4.6 The secondary mode and the frequency spacing 

While the frequency septuplet of the primary mode is consistent 
with our model and values of i and /3, the secondary frequency 
triplet of V2 is problematic. From the data in Table 4 we find for v\ 
that + A'~"1'')/A\^^'^ = 3.96 ± 0.07, which is significantly 



l(»'2) 



1.51 ± 



different to the value for V2 of {A'^+l' + A^^l')/A)^ 
0.23. This indicates that the two modes either arise from different 
spherical degrees, I, or have different pulsation axes. 

We filtered all but the V2 triplet from the data, then calculated 
amplitude and phase modulation of that with respect to rotation 
phase. Because of the smaller amplitude, we needed to use 20 cy- 
cles for each fit, rather than the 5 we used for v\. Fig. 15 shows the 
results. Amplitude maximum occurs at rotation phase 0.67. There is 
an appearance of perhaps another maximum around rotation phase 
0.15; this is probably real, since the rotational sidelobes sum to 
more than the amplitude of the central peak of the triplet, hence 
two maxima per rotation are expected. That the rotational phase of 
amplitude maximum does not coincide with that for supports 



Table 7. This table is the same as the part of Table 4 concerning V2, except 
the zero point in time has been chosen to give phase equality for the com- 
ponents of the U2 multiplet. This time zero point is BJD 245 5172.70203; 
it differs from the time of maximum for vi by 3.7902 d, or 0.67 rotation 
periods. The other frequencies were included in the fit, but are not shown 
here, as only their phases have changed from Table 4; they are no longer 
infonnative with this different time zero point. 
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frequency 


amplitude 


phase 




/xHz 


Mmag 


radians 


V2 - VroX 


917.5126 


8.6 ± 1.3 


3.02 ±0.15 


V2 


919.5464 


11.6 ± 1.3 


3.00 ±0.11 


V2 + VroX 


921.5801 


8.6 ± 1.3 


3.02 ±0.15 




0,0 0.5 1.0 1.5 

rotation phase 



2.0 



Figure 15. Top: This shows the amplitude modulation as a function of rota- 
tional phase for V2 ■ It has been created by fitting 20 cycles of data by linear 
least squares for V2 ■ Notice that the amplitude is modulated only once per 
cycle; both amplitude and phase are poorly defined away from amplitude 
maximum, which occurs around rotation phase 0.67. 



the suggestion that the two modes are either of different spherical 
degree, or have different pulsation axes. We also re-ran the least 
squares fit of Table 5, but shifted the zero point in time to get phase 
equality for the P2 triplet, as shown in Table 7. All three frequen- 
cies show the same phase within the errors for a time of maximum 
at rotational phase 0.67. This is consistent with Fig. 15. 

The separation of vi and i/2 is vi — 1^2 = 55.07 ± 0.03 pHz. 
The model presented in Section 4.2 has a large separation Avq ~ 
57 /iHz, which is close enough to the observed separation to sug- 
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gest that U2 is also a dipole mode. If that is true, then the two modes 
do not have the same pulsation axis. On the other hand, if we dis- 
miss the model and conclude that the degrees of the two modes are 
of different parity, then 55 ^Hz cannot be the large separation. 

From Tcff = 7200 K and i? = 3.6 R©, we have logTefi = 
3.86 and logL/L© = 1.5. From a grid of non-magnetic models 
constructed using the stellar evolution code CESAM (Morel 1997) 
and corresponding pulsation frequencies computed with the code 
ADIPLS (Christensen-Dalsgaard 2008), we find that a model with 
these global properties has a large separation of 25 /iHz. On the 
other hand, a large separation of about 2/3 of the observed value is 
realized by a slightly less luminous model, with log I//Lq = 1.3, 
and similar effective temperature. With the case for vi being a 
dipole mode, a plausible alternative is thus that 1^2 is a radial or 
quadrupole mode and that the star is more luminous than what 
would be expected by assuming the observed separation to be the 
large separation, having, instead, Ai^o ~ 37 /iHz. In this case the 
two observed modes are not consecutive in radial overtone; i.e., 
there is a mode between them that is not excited to observable am- 
plitude. This is seen in another roAp star, HD 217522 (Kreidl et al. 
1991), so is not without precedent. The difficulty with this solution, 
when considering the effect of a dipolar magnetic field on the oscil- 
lations, is that the magnetic distortion of the eigenfunctions of even 
degree modes always generates w I = Q component. Even when 
this component is comparably smaller than the I = 2 component, 
after averaging over the disk its impact is significant. In particu- 
lar, we were not able to find a magnetic model within the scenario 
of V2 being of degree I = 2, 'm which the pulsation phase of this 
second mode jumps by an amount comparable to that seen in the 
observations at particular rotation phases. 

Fig. 15 shows that maximum amplitude for 1^2 is at rotation 
phase 0.67; this is definitely different from the maximum phase of 
v\. Thus the rotational phase of maximum and the multiplet am- 
plitude ratios for u\ and 1/2 are clearly different, yet y\ — V2 = 
55.07 /^Hz is a reasonable value as a separation of the two consec- 
utive overtones, suggesting that the modes are of the same degree, 
I. A first approach is to treat V2 as if it were a pure oblique dipole 
mode perpendicular to //i: 72 ~ 7r/2 + 71. The maximum ampli- 
tude phase (~ 0.67) and the sub-peak amplitude phase (~ 0.15) 
coincide with extrema of the rotational phase, thus the interpreta- 
tion of the V2 mode as an oblique dipole mode is reasonable. In that 
case, the amplitude ratio is 

"^+1 +^-1 ^ tan^ 
^(-2) ~ tan 71- 

However, with the values of the inclinations of the mode Pi and 
the observer (i — ()2^ circ, 71 — 63.5°, derived in section 4.3), 
we cannot fit the amplitude ratio of Eq. 2 within the error bars. We 
conclude that the dipole mode V2 is neither aligned nor perpendic- 
ular to i/i . In Section 5 a tentative explanation is proposed which 
suggests that the magnetic field acts in a different way on the two 
modes. 



4.7 Asteroseismic position in the HR Diagram 

The position of KIC 10195926 in the Hertzsprung-Russell (HR) 
Diagram can be compared with the theoretical instability strip for 
roAp stars modelled by Cunha (2002). The computation of the in- 
stability strip, shown in Fig. 16, was based on a pulsational stability 
analysis of models which details are described in Balmforth et al. 
(2001). These are non-magnetic models, in which the effect of the 



magnetic field enters only in an indirect way, through an assumed 
suppression of envelope convection. 

The stability against high radial order pulsations of a number 
of main-sequence models of different masses was considered. For 
each model with unstable high radial order modes, the frequency 
of the mode with largest growth rate was taken as a reference for 
the typical frequencies that may be expected to be observed in that 
region of the HR diagram. The reference frequencies are shown in 
Fig. 16 where we overplotted the three models discussed in the pre- 
vious sections. It is clear that the reference frequencies decrease as 
the star evolves, a phenomenon that is in great part (but not totally) 
explained by the increase in stellar radius (see Cunha 2002, for a 
detailed discussion). 

From Fig. 16 it can be seen that the star is near the terminal age 
main sequence, as expected from its low effective temperature and 
surface gravity. That provides a natural explanation for the low fre- 
quency of the modes observed in KIC I0I95926, when compared 
to those seen in the majority of roAp stars. 

The models considered in the previous sections, as well as 
those used to produce the instability strip shown in Fig. 16, all have 
metallicity Z = 0.02. On the other hand, our analysis of the spectra 
of KIC 10195926 in Section 2 indicates that its photospheric mean 
iron abundance is comparable to or slightly higher than the corre- 
sponding solar abundance of Asplund et al. (2009). As we pointed 
out, the global abundances may be greater than solar, given the 
star's young age compared to the Sun and galactic chemical evo- 
lution. Thus our adopted metallicity in our models of Z = 0.02 
seems a reasonable first choice. 

The effects of metallicity on the models need to be explored. 
This will be done when more detailed spectroscopic analysis of 
the star becomes available. Nevertheless, we can anticipate that the 
main impact of a variation in the metallicity adopted in the models 
will be that of a change of the mass derived for the best model of 
KIC 10195926. The effect of the metallicity was considered, e.g., in 
the modelling of the roAp starHR 12 17 by Cunha et al. (2003), who 
found that for a fixed stellar radius, changing from Z — 0.009 to 
Z = 0.019 led to a change of about 4 /iHz in Av (which is 68 /iHz 
for HR 1217), as a consequence of a change of ~O.2M0 in the 
stellar mass. These variations, which are also expected if models 
of different metallicity are considered for KIC 10195926, do not 
change any of the conclusions presented in this work. We also note 
that we expect that decreasing the metallicity would hardly change 
the theoretical results presented in Figs. 11 to 13, although, again, 
the parameters of the best model might be slightly modified. 



5 DISCUSSION 

5.1 The pulsation axes of the two modes 

KIC 10195926 shows the power of the precision of the Kepler pho- 
tometric data. While its principal pulsation at Pi could have been 
discovered in ground-based observations, only the /imag precision 
of the data has allowed us to see examine the rotational multiplet 
of i/i, with all the mode geometry information that contains. And 
z^2 could not be detected from ground-based photometry. It is the 
separation between, and comparison of, Pi and P2 that lead to the 
remarkable suggestion that the pulsation axes of these two modes 
do not coincide. This is a viable hypothesis for this star supported 
by the improved oblique pulsator model (Bigot & Dziembowski 
2002). No other pulsating star has ever been noted to have different 
pulsation axes for different modes. 
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Figure 16. HR diagram with the position of the models discussed in Sections 4.2 and 4.6. The full lines show evolutionary tracks for masses between 1 .5 Mq 
and 2.4 Mq. The numbers along the tracks show the frequency, in mHz, of the most unstable mode for each model considered (marked by open circles). 
The large filled circles show models that are stable against pulsations for frequencies in the roAp star range (adapted from Cunha 2002). Dashed lines show 
the ranges of effective temperature (vertical) and luminosity (oblique) expected for the star from the non-seismic observations alone. The luminosity limits 
were derived from the grid of non-magnetic models, assuming 3.5 ^5 log g ^ 4.0. The positions of three stellar models are shown in this diagram, namely, 
the magnetic model discussed in Section 4.2 (cross), the model discussed in Sections 4.4 and 4.6 (diamond) derived from non- seismic data alone, and the 
non-magnetic model with large separation of ^ 37/iHz, discussed in Section 4.6 (star). 



In the simple oblique pulsator model (Kurtz 1982) the pulsa- 
tion axis of the star coincides with its magnetic axis. The simple 
case for the rotational light variations of CVn stars is also for 
spots that are concentric about the magnetic and pulsation poles. It 
has been clear for decades that magnetic Ap stars often have global 
fields that are more complex than dipoles, and that the spot struc- 
ture is more detailed than two spots at the poles. Recently, Bruntt 
et al. (2009) showed that the spots of the roAp star a Cir do not co- 
incide with its pulsation axis. We now see for KIC 10195926 that 
u\ rotational light minimum in brightness coincides with maximum 
pulsation amplitude, whereas for V2 maximum pulsation amplitude 
coincides with one of the rotational maxima in brightness. These 
indicate that the modes are in alignment with spots, but not with 
the same ones, and there is no simple symmetry to the spot struc- 
ture. 

No other roAp star so obviously shows these effects. The num- 
ber of stars with rotational multiplets is small. HD6532, HR3831 
and HR 1217 are the best cases. For the first two there is only 
one pulsation mode, while the latter shows at least 6 modes that 
have been modelled with distorted modes of different spherical 
degree, but with no thought to use differing pulsation axes (Saio 
et al. 2010). In light of the result here for KIC 10195926, it will be 



worth revisiting the interpretation of the frequency multiplets for 
HR1217. 

The theory of the interaction of rotation, pulsation and mag- 
netic fields is complex. Bigot & Dziembowski (2002) presented the 
improved oblique pulsator model incorporating both rotation and 
magnetic field in the relatively weak field case where Bp ~ 1 kG. 
They showed that the pulsation axis is not necessarily the magnetic 
axis of the star. In their formalism, the inclination of the dipole 
mode depends on how the centrifugal frequency shift compares to 
the difference between magnetic eigenmode frequencies of consec- 
utive azimuthal orders m, i.e. A;;'f = uj'^'^if^ - 

At a given rotation rate, the difference between magnetic and 
centrifugal shifts is a function of the magnetic strength and of 
the frequency. Therefore, strictly speaking two consecutive mag- 
netic overtones of the same degree / must have different inclina- 
tions with respect to the magnetic field axis. However, as shown in 
Bigot et al. (2000), the difference A™^^ — ^"±1 i ^^"^y small in 
most cases, which indicates that consecutive modes have roughly 
the same inclination. Since v\ and V2 have different amplitude ra- 
tios we conclude that the two modes have different axes. One pos- 
sibility to explain this is to consider that the two modes are perpen- 
dicular. 

Since we found that the vi mode is a dipole almost linearly po- 
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larized close to the magnetic axis, a simple interpretation for hiav- 
ing a mode at V2 with a different pulsation axis would have been 
an orthogonal dipole mode. However, we could not find satisfac- 
tory agreement in fitting both amplitude ratios simultaneously for 
u\ and U2, which shows that this hypothesis has to be discarded. A 
mode of degree Z = 3 is also a possibility, but the phase reversal 
of V2 occurring between times of pulsation amplitude maxima, as 
seen in Fig. 15, and as deduced from the amplitudes of the V2 fre- 
quency triplet given in Table 7, suggests the simple geometry of a 
dipole. 

Another possible explanation for the difference in the pulsa- 
tion axes of the two modes observed could come from an effect 
found by Cunha & Gough (2000) and Saio & Gautschy (2004). 
They clearly showed that at specific frequencies, which depend on 
the magnetic field strength, two modes of consecutive radial order 
are very differently affected by the magnetic field. In those cases, 
the relative size of magnetic and centrifugal shifts is likely to be 
strongly modified from one overtone to the next, and the pulsation 
axes of the two consecutive modes could be very different. Further 
investigation needs to be carried out to test this possibility. 



5.2 The rotational subharmonic 

An important issue is the presence of the rotational subharmonic 
which manifests itself as a light variation with twice the rotational 
period; see Fig. 5. This is clearly shown by the presence of the sub- 
harmonic of the orbital frequency at \ in the Fourier Transform 
of independent light curves. We could take the view that this varia- 
tion is external to the star. We could imagine, for example, a com- 
panion in an orbit with twice the rotational period and exerting a 
2: 1 resonant tidal effect, in which case none of the conclusions so 
far discussed are affected. However, this seems improbable. It will 
be tested with radial velocity measurements in spectra that we are 
obtaining. 

Let us look for a possible cause of the rotational subharmon- 
ics. In what follows we note that the theoretical subharmonics ap- 
pear at = Vrot/'i, 2Vrot/^, Vrotl% 2Vrot/^, iv^ot / ■ ■ ; aS 

well as 4!/iot/3, St'rot/S and 2v^oi- These are a systematic se- 
quence. The observed subharmonics are clearly seen at i^rot /2 and 
St'rot/S. However, least-squares fitting of other possible subhar- 
monics shows that some of them have formally significant ampli- 
tudes. Because of the low frequency noise in the data caused by 
instrumental drift - particularly in Ql - and because a gap in the 
data between Ql and Q3.3 introduces aliases, we are not confident 
in the formal errors. In what follows the frequencies proposed are 
consistent with the data, with v-^ot /2 and S/^rot /2 having the high- 
est confidence. 

Our hypothesis is that torsional modes, or r modes, are excited 
and generate the observed subharmonic frequencies. In a rotating 
star, the radial component of vorticity interacts with the Coriolis 
force and disturbs the equilibrium structure of the star This gener- 
ates r-mode oscillations, which are torsional oscillations (Unno et 
al. 1989). 

As the light curve varies with the rotation period, there are 
brightness inhomogeneities on the stellar surface, which may be 
decomposed into a series of spherical harmonics. The mean flux 
distribution averaged over latitude is then expressed in the corotat- 
ing frame as 



where a„i/ denotes the amplitude of brightness inhomogeneity of 
m'- component, and i^r is the azimuthal angle in the corotating 
frame. 

By treating such a patchy configuration as the equilibrium 
state and considering perturbations on it, we expect the flux vari- 
ation caused by the r mode with the angular degree I and the az- 
imuthal order m is expressed as 



5F_ 

F 



(4) 



where Vr denotes the frequency of the r mode with the degree I and 
the azimuthal order m in the corotating frame and it is given, to the 
first order, as 



2mi/rot 

i{i + iy 



(5) 
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F oc 



1 + ^ a„ 



• exp(im'(j!>R) 



(3) 



If m = Z, i/r is reduced to 2z^rot /{I + 1): 2i^rot/3 for I — 2; !/rot 
for/ = 3;2j/rot/5forZ = 4, ... 

In order to move from the corotating frame to an inertial 
frame, we only have to utilize the following transformation; = 
— Vroxt, where 4>i is the azimuthal angle in the inertial frame. 
The observable luminosity variation is then given by 

5L oc exp[im(j!)i -f i2it{v^ — mi^rot)^} 
+ ^ [am' exp{i(m + m')0i 

-\-i2-n{vr - (m + m')frot}t]] ■ (6) 

The component of m' = —m leads the corotating-frame frequency 
of the r mode, i^r, to be observable. The visibility of the corotating- 
frame frequencies depends on the size and the contrast of the sur- 
face brightness inhomogeneity. 

Note that both the corotating-frame and the inertial-frame fre- 
quencies are observable. As an example, we fitted by least-squares 
4urot/S = 0.234552 d"^ to the Q3.3 data and found that it has 
a highly significant amplitude of 52 ± 2 ^mag. In our interpre- 
tation this is a manifestation of the inertial-frame frequencies, 
mi^iot [1 — 2/1(1 + 1)] for I = m = 2. These considerations suggest 
that our new finding of the subharmonics is essentially discovery 
of torsional modes. New data are being acquired by Kepler for this 
star, with which we will test this idea further. 



6 FUTURE WORK 

Kepler is continuing to observe KIC 10195926 in short cadence, 
so we will have a longer data set with higher S/N to work with 
in the future. We have obtained high resolution spectra and will 
gather more in future observing seasons. These will provide better 
constraints on T^s and logg, which because of their current er- 
rors leave uncertainties in the interpretation of and 1/2 within our 
models, hence uncertainty in the large separation. KIC 10195926 
is clearly a spectrum variable, and we plan Doppler Imaging stud- 
ies of its spectral line variations to map surface abundance distribu- 
tions; we also plan model reconstruction of the spots that lead to the 
rotational photometric variations. Spectropolarimetry is planned to 
determine the magnetic field strength and structure. KIC 10195926 
is an outstanding example of the power of the time span, duty cy- 
cle and precision of the Kepler data for asteroseismology, and it 
demonstrates well the need and usefulness of ground-based follow- 
up observations using many medium aperture telescopes. 
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